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Abstract 

Prognostics,  which  deals  with  predicting  remaining  useful 
life  of  components,  subsystems,  and  systems,  is  a  key  tech¬ 
nology  for  systems  health  management  that  leads  to  improved 
safety  and  reliability  with  reduced  costs.  The  prognostics 
problem  is  often  approached  from  a  component-centric  view. 
However,  in  most  cases,  it  is  not  specifically  component  life¬ 
times  that  are  important,  but,  rather,  the  lifetimes  of  the  sys¬ 
tems  in  which  these  components  reside.  The  system-level 
prognostics  problem  can  be  quite  difficult  due  to  the  increased 
scale  and  scope  of  the  prognostics  problem  and  the  rela¬ 
tive  lack  of  scalability  and  efficiency  of  typical  prognostics 
approaches.  In  order  to  address  these  issues,  we  develop 
a  distributed  solution  to  the  system-level  prognostics  prob¬ 
lem,  based  on  the  concept  of  structural  model  decomposi¬ 
tion.  The  system  model  is  decomposed  into  independent 
submodels.  Independent  local  prognostics  subproblems  are 
then  formed  based  on  these  local  submodels,  resulting  in  a 
scalable,  efficient,  and  flexible  distributed  approach  to  the 
system-level  prognostics  problem.  We  provide  a  formulation 
of  the  system-level  prognostics  problem  and  demonstrate  the 
approach  on  a  four-wheeled  rover  simulation  testbed.  The  re¬ 
sults  show  that  the  system-level  prognostics  problem  can  be 
accurately  and  efficiently  solved  in  a  distributed  fashion. 

1.  Introduction 

Prognostics  is  the  process  of  predicting  the  end  of  (useful)  life 
(EOL)  and/or  the  remaining  useful  life  (RUL)  of  components, 
subsystems,  or  systems.  The  prognostics  problem  itself  can 
be  divided  into  two  distinct  problems:  (i)  the  estimation  prob¬ 
lem,  which  determines  the  current  state  of  the  system,  and  (if) 
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the  prediction  problem,  which,  using  the  current  system  state 
estimate,  computes  EOL  and/or  RUL.  In  this  paper,  we  focus 
on  a  model-based  prognostics  approach  (Orchard  &  Vachtse- 
vanos,  2009;  Daigle  &  Goebel,  201  lb;  Saha  &  Goebel,  2009; 
Luo  et  al.,  2008).  In  model-based  prognostics,  an  underly¬ 
ing  model  of  the  system,  its  components,  and  how  they  fail 
is  leveraged,  where  health  state  estimation  is  formulated  as 
a  joint  state-parameter  estimation  problem,  typically  using  a 
filtering  approach,  and  prediction  is  formulated  as  a  simula¬ 
tion  problem  (Daigle,  Saha,  &  Goebel,  2012). 

To  the  best  of  our  knowledge,  all  prognostics  research  to  date 
has  been  focused  on  individual  components,  and  determining 
their  EOL  and  RUL,  e.g.,  (Orchard  &  Vachtsevanos,  2009; 
Saha  &  Goebel,  2009;  Daigle  &  Goebel,  2011a;  Celaya  et 
al.,  2011;  Bolander  et  al.,  2010;  Luo  et  al.,  2008;  Bying- 
ton  et  al.,  2004).  However,  in  many  cases,  the  desired  infor¬ 
mation  is  the  EOL  of  the  system,  which  is  obtained  through 
system-level  prognostics.  Generally,  the  EOL  of  a  system  de¬ 
pends  on  its  constituent  components  and  how  they  interact. 
Approaching  this  problem  from  the  centralized  perspective 
becomes  very  difficult,  as  common  (centralized)  prognostics 
algorithms  may  not  scale  to  the  system  level. 

In  order  to  address  the  problems  with  centralized  approaches, 
in  recent  work,  we  have  developed  a  distributed  model-based 
prognostics  architecture  that  allows  the  decomposition  of  a 
large  prognostics  problem  into  several  independent  local  sub¬ 
problems  from  which  local  results  can  be  merged  into  a  global 
result  (Daigle  et  al.,  2011;  Daigle,  Bregon,  &  Roychoudhury, 
2012).  Since  each  local  subproblem  can  be  solved  indepen¬ 
dently,  each  can  be  assigned  to  a  different  processing  unit 
and  be  solved  in  parallel.  Such  a  distributed  approach  is  in 
contrast  to  other  proposed  distributed  prognostics  architec¬ 
tures  in  which  the  global  problem  is  not  decomposed  and 
the  computation  is  distributed  onto  multiple  processing  units, 
e.g.,  (Saha,  Saha,  &  Goebel,  2009).  Our  distributed  approach 
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scales  well  and  the  resulting  subproblems  are  typically  small 
and  easy  to  solve,  resulting  in  an  efficient  and  flexible  dis¬ 
tributed  solution  to  the  prognostics  problem.  Such  an  ap¬ 
proach  has  obvious  advantages  when  applied  to  the  system- 
level  prognostics  problem.  In  this  paper,  we  formulate  the 
system-level  prognostics  problem  and  propose  a  solution  us¬ 
ing  this  distributed  prognostics  framework.  We  apply  our 
system-level  prognostics  approach  to  a  rover  testbed  and  pro¬ 
vide  results  in  simulation  to  empirically  demonstrate  and  val¬ 
idate  the  approach. 

The  paper  is  organized  as  follows.  Section  2  formulates 
the  system-level  prognostics  problem  and  overviews  the  pro¬ 
posed  distributed  solution.  Section  3  describes  the  estima¬ 
tion  problem,  and  Section  4  describes  the  prediction  problem. 
Section  5  presents  the  rover  case  study,  and  shows  prognos¬ 
tics  results  in  simulation.  Section  6  concludes  the  paper. 

2.  System-level  Prognostics 

While  most  prognostics  approaches  focus  on  individual  com¬ 
ponents,  in  most  practical  cases  it  is  actually  the  EOL  of  the 
system  that  must  be  determined.  With  this  prediction,  the 
future  usage  of  the  system  may  be  optimally  planned  to  max¬ 
imize  system  life  and  to  schedule  system-wide  maintenance 
activities.  It  is  often  important  to  take  a  system-level  per¬ 
spective  of  prognostics,  because  the  degradation  of  individual 
components  is  often  coupled,  i.e.,  the  way  one  component  de¬ 
grades  is  dependent  on  how  a  connected  component  degrades. 
This  may  occur,  for  example,  if  one  component  provides  the 
inputs  to  another  component,  in  which  case,  prognostics  of 
the  latter  component  cannot  be  performed  in  isolation. 

In  this  section,  we  first  define  the  system-level  prognostics 
problem.  We  then  introduce  the  system-level  prognostics  ap¬ 
proach  and  architecture  using  a  distributed  prognostics  frame¬ 
work  that  is  based  on  structural  model  decomposition. 

2.1.  Problem  Formulation 

The  goal  of  system-level  prognostics  is  the  prediction  of  the 
EOL  and/or  RUL  of  a  system.  We  assume  the  system  model 
may  be  generally  defined  as 

x(i)  =  f(f,x(f),0(f),u(f),v(f)), 
y (t)  =  h(t,  x(t),  9(f),  u(t),  n(t)), 

where  x(t)  £  R"x  is  the  state  vector,  9(f)  £  R”9  is  the 
unknown  parameter  vector,  u (t)  £  R"“  is  the  input  vector, 
v(t)  £  M””  is  the  process  noise  vector,  f  is  the  state  equation, 
y(t)  £  R"9  is  the  output  vector,  n(t)  £  K""  is  the  measure¬ 
ment  noise  vector,  and  h  is  the  output  equation. 1  This  model 
describes  both  the  nominal  behavior  and  faulty  behavior,  in¬ 
cluding  the  fault  progression  functions. 


'Here,  we  use  bold  typeface  to  denote  vectors,  and  use  na  to  denote  the 
length  of  a  vector  a. 


In  system-level  prognostics,  we  are  interested  in  when  the 
performance  of  a  system  lies  outside  some  desired  region 
of  acceptable  behavior.  The  desired  performance  is  ex¬ 
pressed  through  a  set  of  nc  constraints,  Ceol  =  {cj}”=  i> 
where  c,  :  Rnx  x  R™9  x  Rn“  — y  B  maps  a  given  point 
in  the  joint  state-parameter  space  given  the  current  inputs, 
(x(f),  9(f),  u (/)),  to  the  Boolean  domain  B  =  [0, 1],  where 
Cj(x(f),  0(f),  u(f))  =  1  if  the  constraint  is  satisfied.  If 
Ci(x(t),  9(t),u(t))  =  0,  then  the  constraint  is  not  satis¬ 
fied,  and  the  behavior  of  the  system  is  deemed  to  be  un¬ 
acceptable.  These  deterministic  constraints  may  refer  to 
component-level,  subsystem-level,  or  system-level  specifica¬ 
tions  or  requirements  and  define  a  fixed  partition  of  the  state- 
parameter-input  space  into  acceptable  and  unacceptable  re¬ 
gions  of  behavior.  When  the  constraints  are  violated,  it  does 
not  necessarily  refer  to  a  hard  failure,  but  any  point  at  which 
the  operational  risk  is  too  large  to  continue  system  operation, 
or  future  behaviors  of  the  system  will  be  in  some  way  unac¬ 
ceptable.  At  this  point  we  say  the  system  has  no  useful  life 
remaining. 

These  individual  constraints  may  be  combined  into  a  single 
system-level  threshold  function  Teol  ■  Rrax  x  R"9  x  Rn“  — > 
B,  defined  as 

TEOL(x-{i),9(t),  u(f))  = 

1»  0  g  {c;(x(f),0(f),u(*))}"=i 

0,  otherwise. 

Teol  evaluates  to  1,  i.e.,  the  system  has  reached  an  unac¬ 
ceptable  region  of  behavior,  when  any  of  the  constraints  are 
violated.  EOL  is  then  defined  as 

EOL[tP)  = 

inf{i  g  R  :  t  >  tP  A  TEOL(x(f),9(f),  u(f))  =  1}, 

i.e.,  EOL  is  the  earliest  time  point  at  which  Teol  is  met  (eval¬ 
uates  to  1).  RUL  is  expressed  using  EOL  as 

RUL{tP)  =  EOL(tP)  —  tP. 

Note  that  because  x(f)  is  a  random  variable,  EOL  and  RUL 
must  necessarily  be  random  variables  also. 

2.2.  Prognostics  Approach 

In  order  to  make  an  EOL  or  RUL  prediction  for  the  system, 
the  initial  state  from  which  to  make  a  prediction  is  required. 
In  general,  this  initial  state  is  not  directly  observed,  and  must 
be  estimated.  Therefore,  there  are  two  sequential  problems 
for  prognostics:  the  estimation  problem  and  the  prediction 
problem.  The  estimation  problem  is  to  find  a  joint  state- 
parameter  estimate  p(x(t),  9(t)\y0:t)  based  on  the  history  of 
observations  up  to  time  t,  yo-.t-  This  estimate  is  represented 
as  a  probability  distribution  because,  generally,  the  system 
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state  is  not  directly  observed,  and  there  is  sensor  noise,  n(£), 
and  process  noise,  v(f).  At  a  given  prediction  time,  tp ,  the 
prediction  algorithm  uses  the  joint  state-parameter  estimate 
P(x(tp),0(tp) |y0:tP)  and  computes p(EOL(tP)\y0-.tP)  and 
p(RU L(tp)\yo-.tp)-  Along  with  the  uncertainty  in  the  state- 
parameter  estimate,  process  noise  and  uncertainty  in  the  fu¬ 
ture  inputs  to  the  system  all  contribute  to  the  uncertainty  in 
the  EOL/RUL  prediction. 

This  system-level  prognostics  problem,  consisting  of  estimat¬ 
ing  the  system  state  and  then  predicting  its  evolution  to  EOL, 
can  be  solved  using  component-level  approaches  by  treating 
the  entire  system  as  a  single  component  and  applying  these 
approaches  directly.  However,  for  a  large  system,  both  the  es¬ 
timation  and  prediction  problems  are  correspondingly  large. 
Due  to  the  large  state -parameter  dimension,  a  centralized  ap¬ 
proach  does  not  scale  well,  and  can  be  very  inefficient. 

Therefore,  we  propose  to  decompose  the  global  system- 
level  prognostics  problem  into  independent  local  subprob¬ 
lems,  such  that  the  solutions  to  the  local  subproblems  may 
be  easily  merged  to  form  the  solution  to  the  global  prognos¬ 
tics  problem.  This  forms  a  naturally  distributed  approach  in 
which  the  local  subproblems,  since  they  are  independent,  may 
be  solved  in  parallel,  thus  providing  scalability  and  efficiency. 
Further,  the  approach  allows  different  algorithms  to  be  em¬ 
ployed  on  each  subproblem.  The  subproblems  often  corre¬ 
spond  directly  to  component-level  prognostics  problems,  and 
the  approach  provides  a  mechanism  to  combine  component- 
level  prognostics  results  into  system-level  results. 

In  (Daigle  et  al.,  201 1),  we  developed  such  a  distributed  solu¬ 
tion  to  the  estimation  part  of  the  prognostics  problem,  based 
on  the  concept  of  structural  model  decomposition  (Pulido  & 
Alonso-Gonzalez,  2004).  In  recent  work,  the  same  concept 
was  used  to  decompose  the  prediction  problem  (Daigle, 
Bregon,  &  Roychoudhury,  2012).  Structural  model  decom¬ 
position  allows  one  to  decompose  a  system  model  into  a  set 
of  submodels  for  which  local  prognostics  problems  can  be 
directly  defined.  The  global  model  of  the  system,  denoted  as 
Ad,  is  defined  as  follows. 

Definition  1  (Model).  The  model  of  a  system.  Ad,  is  a  tuple 
Ad  =  (X,  ®,U,Y,C ),  where  X  is  the  set  of  state  variables 
of  x,  0  is  the  set  of  unknown  parameters  of  0,  U  is  the  set  of 
input  variables  of  u,  Y  is  the  set  of  output  variables  of  y,  and 
C  is  the  set  of  model  constraints  of  f,  h,  and  Cpo  l  ■ 

Informally,  a  model  consists  of  a  set  of  variables  and  a  set 
of  constraints  among  the  variables.  While  technically  f  and 
h  themselves  are  (complex)  constraints,  we  represent  them 
instead  as  sets  of  simple  constraints.  This  view  is  also  more 
consistent  with  the  way  modelers  describe  f  and  h,  i.e.,  as  sets 
of  equations,  each  describing  a  single  state  or  output  variable. 

Model  decomposition  is  accomplished  by  assigning  some 
variables  as  local  inputs  for  which  the  values  are  known  (e.g.. 


they  are  directly  measured).  In  this  way,  the  submodels  are 
made  computationally  independent  of  each  other.  Within  this 
scheme,  a  submodel  is  then  defined  as  follows. 

Definition  2  (Submodel).  A  submodel  Adi  of  a  system  model 
Ad  =  (X,  0,  U,  y,  C)  is  a  tuple  M,  =  {Xt,  U„  Yu  C'  ), 
where  C  X,  0,  C  0,  U,  C  X  U  U  U  Y,  and  F,  C  T  are 
the  state,  parameter,  input,  and  output  variables,  respectively, 
and  Ci  C  C  are  the  submodel  constraints. 

For  distributed  prognostics,  we  find  a  set  of  submodels  that 
satisfy  a  certain  set  of  properties.  For  distributed  estimation, 
the  submodels  use  Ul  C  U  U  [Y  —  Yj),  and  we  find  a  set  of 
minimal  submodels  such  that  each  Tj  is  a  singleton,  and  over 
all  Yi,Yj  where  i  ^  j,  Yi:  fi  Y)  =  0.  So,  each  submodel  uses 
some  global  model  inputs  and  some  measured  values  as  lo¬ 
cal  inputs,  and,  in  this  way,  the  submodels  become  decoupled 
and  may  be  computed  independently  from  each  other.  By  cre¬ 
ating  submodels  with  one  output  variable  each,  we  maximize 
the  number  of  estimation  submodels  and  the  opportunity  for 
parallelization  of  the  estimation  task.  By  making  the  sub¬ 
models  minimal,  they  require  no  constraints  or  variables  that 
are  not  strictly  necessary  to  compute  Y,.  An  algorithm  for 
computing  the  set  of  submodels  with  these  properties  is  given 
in  (Daigle  et  al.,  2011),  which  is  based  on  the  model  decom¬ 
position  algorithms  presented  in  (Pulido  &  Alonso-Gonzalez, 
2004;  Bregon,  Biswas,  &  Pulido,  2012). 

For  distributed  prediction,  the  submodels  use  Ut  C  Up, 
where  Up  C  X  U  U.  Here,  Up  is  a  set  of  variables  whose 
future  values  can  be  hypothesized.  In  the  centralized  case. 
Up  =  U.  We  find  a  set  of  minimal  submodels  such  that 
each  submodel  has  at  least  one  c  £  Ceol  belonging  to  Ci, 
and  over  all  submodels,  Ceol  is  covered.  This  ensures  that 
Teol  may  be  computed  for  the  system;  since  Teol  is  1 
whenever  any  of  the  constraints  in  Ceol  are  violated,  we  can 
independently  evaluate  when  those  individual  constraints  will 
be  violated  and  then  take  the  minimum  to  obtain  the  system 
EOL.  An  algorithm  for  computing  the  set  of  submodels  with 
these  properties  is  given  in  (Daigle,  Bregon,  &  Roychoud¬ 
hury,  2012).  Both  decomposition  algorithms  work  in  a  sim¬ 
ilar  way;  essentially,  they  start  with  a  variable  or  constraint 
that  must  be  computed  in  the  local  submodel,  and  then  trace 
the  dependencies  backwards  until  local  inputs  are  reached, 
including  all  variables  and  constraints  found  throughout  the 
search  within  the  submodel. 

Note  that  the  problem  of  defining  Up  is  critical  to  obtaining 
accurate  results  for  system-level  EOL  in  a  distributed  manner. 
On  average,  the  most  accurate  result  will  be  achieved  when 
the  system  model  is  directly  used  for  prediction,  because  it 
captures  all  the  interdependencies  between  the  components. 
In  the  general  case,  damage  could  be  progressing  in  multiple 
components,  and  damage  progression  in  one  component  may 
have  an  effect  on  damage  progression  in  another  component 
due  to  their  coupling.  In  such  cases,  for  system-level  prog- 
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Prognostics 


Figure  1 .  Sample  system-level  prognostics  architecture. 


nostics  the  components  cannot  be  decoupled  due  to  these  in¬ 
teractions,  and  the  prediction  problem  cannot  be  decomposed 
into  two  independent  problems,  one  for  each  component.  It  is 
only  appropriate  to  neglect  these  interactions  when  they  are 
either  negligible  or  predictable  a  priori.  It  will  be  shown  in 
Section  5  how  this  is  an  important  consideration. 

2.3.  Prognostics  Architecture 

A  sample  system-level  prognostics  architecture  based  on  the 
distributed  framework  is  shown  in  Fig.  1.  In  discrete  time 
k ,  and  using  a  discrete-time  version  of  the  model,  the  dam¬ 
age  estimation  module  takes  as  input  both  uk  and  y/,:  and 
splits  them  into  local  inputs  and  outputs  for  the  submodels. 
Estimation  is  performed  for  each  submodel  using  an  appro¬ 
priate  algorithm,  computing  local  state-parameter  estimates 
p(xj.,  Olk\yQ.k).  Some  of  these  local  estimates  are  merged 
corresponding  to  the  prediction  submodels.  For  example, 
submodel  M5  builds  its  local  state  using  the  estimates  from 
the  estimators  of  M2,  M3,  and  M 4.  The  local  predictors 
compute  local  EOL/RUL  predictions  p(EOLkp\yl0.kp)  and 
p(RULkp\yl0.kp)  at  given  prediction  time  kp  based  on  the 
local  EOL  constraints.  Local  predictions  are  then  merged  into 
global  predictions  p(EOLkp\y0:kp)  and  p(RULkp  |y0:fcp) 
by  taking  the  minimum  of  the  local  predictions. 

3.  Distributed  Estimation 

As  described  in  Section  2,  in  our  distributed  estimation 
scheme,  the  local  estimator  for  each  submodel  Mi  produces 
a  local  estimate  p(yLlk,Olk |yp.fc),  where  xj,  C  xfc,  0k  C  9k, 
and  y\  C  yk.  Here,  the  local  inputs  used,  u\  consist  of  ele¬ 
ments  from  both  u  and  y,  where  measured  values  are  directly 
used  as  local  inputs.  The  estimation  problem  is  decomposed 
by  finding  a  set  of  minimal  submodels  that  together  cover  the 
subset  of  x  and  6  required  for  prediction,  by  using  these  local 
inputs.  This  approach  to  distributed  estimation  is  different 
from  approaches  like  the  distributed  decentralized  extended 
Kalman  filter  (Mutambara,  1998)  or  other  estimation  fusion 


techniques  (Sinha  et  al.,  2008)  where  local  estimates  are  com¬ 
municated  between  local  estimators.  Here,  local  estimators 
do  not  communicate  and  operate  completely  independently. 

In  order  to  effectively  perform  joint  state-parameter  estima¬ 
tion,  the  system  should  be  observable,  among  other  require¬ 
ments.  If  the  global  model  is  structurally  observable,  then  we 
are  guaranteed  that  the  local  submodels  for  estimation  are  as 
well  (Moya  et  al.,  2010). 

Any  suitable  algorithm  may  be  used  for  joint  state-parameter 
estimation.  In  this  paper,  we  use  an  unscented  Kalman  fil¬ 
ter  (UKF)  (Julier  &  Uhlmann,  1997,  2004)  with  a  variance 
control  algorithm  (Daigle,  Saha,  &  Goebel,  2012).  The  UKF 
assumes  the  general  nonlinear  form  of  the  state  and  output 
equations  described  in  Section  2,  but  restricted  to  additive 
Gaussian  noise. 

We  summarize  the  main  details  of  the  UKF  below,  and  refer 
the  reader  to  (Julier  &  Uhlmann,  1997,  2004)  for  details.  In 
the  UKF,  distributions  are  approximated  using  the  unscented 
transform  (UT).  The  UT  takes  a  random  variable  x  £  R”x, 
with  mean  x  and  covariance  Pxx,  that  is  related  to  a  second 
random  variable  y  £  by  some  function  y  =  g(x),  and 
computes  the  mean  y  and  covariance  Pyy  using  a  minimal  set 
of  deterministically  selected  weighted  samples,  called  sigma 
points  (Julier  &  Uhlmann,  1997).  X'  denotes  the  /th  sigma 
point  from  x  and  wl  denotes  its  weight.2  The  sigma  points 
are  always  chosen  such  that  the  mean  and  covariance  match 
those  of  the  original  distribution,  x  and  Pxx.  Each  sigma 
point  is  passed  through  g  to  obtain  new  sigma  points  y,  i.e., 

y  =  g(xi) 


2 Sigma  point  weights  do  not  directly  represent  discrete  probabilities,  so  are 
not  restricted  to  [0, 1]. 
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with  mean  and  covariance  calculated  as 


step  is: 


y =Y,wiyl 

i 

P  vv  =  Y,wi{yl-y)(yl-y)T- 


In  this  paper,  we  use  the  symmetric  unscented  transform,  in 
which  2 nx  + 1  sigma  points  are  symmetrically  selected  about 
the  mean  according  to  (Julier  &  Uhlmann,  2004): 


w 


l 


X1 


K 

(' nx  +  k)  ' 

2  (nx  +  k)  ’ 


i  =  0 

i  =  1,.  ■  .,2nx 


p w  —  P-  +  ^2, wl0?k\k-\  ~  yfc|fe-i)(K-|fc-i —  yk\k-if 

i 

ns 

P xy  —  'y  '  w  [X fe|fc_i  —  ^fcl fc —  1 ) (3^fc| fc  —  i  —  y fc| fc  —  1 ) 


i,’  p  r> 

iVfc  r  xy17  yy 

Xfc|fc  =  xfc|fc_!  +  Kfc(yfc  —  yfe|fc— l) 

Pfc|fc  Pfc|fc—  1  PtcPyyPfc  5 

where  R  is  the  sensor  noise  covariance  matrix. 


Joint  state-parameter  estimation  is  accomplished  in  the  UKF 
by  augmenting  the  state  vector  with  the  unknown  parameters, 
and  the  corresponding  diagonal  elements  of  the  process  noise 
matrix,  Q,  are  set  to  nonzero  values.  In  this  way,  the  param¬ 
eter  estimates  become  time-varying  and  are  modified  by  the 
filter  using  the  measured  outputs. 


where  ^  (nx  +  K)Px:2j  refers  to  the  ith  column  of  the  ma¬ 
trix  square  root  of  (nx  +  k) Pxx-  Here,  n  is  a  free  parameter 
that  can  be  used  to  tune  higher  order  moments  of  the  distribu¬ 
tion.  If  x  is  assumed  Gaussian,  then  selecting  k  =  3  —  nx  is 
recommended  (Julier  &  Uhlmann,  1997). 

In  the  filter,  first,  na  sigma  points  Xk-i\k-\  are  derived 
from  the  current  mean  xfc_1|fc_1  and  covariance  estimates 
Pfe_i|fe_i  using  a  sigma  point  selection  algorithm.  The  pre¬ 
diction  step  is: 

^ k\k—l  =  f  (^fc— l|fc— 1’  ufc-l))  i  =  1>  •  ■  •  >  ns 

yk\k-i  =  h(^fc|fc-i)>*  =  1, . . .  ,ns 

ns 

*-k\k—  1  —  ^  ^  W  X h\k—  1 

i 

ns 

Yk\k~i  =  ^2^lylk\k-i 

i 

with 

Pfc|fc-i  =  Q+ 

ns 

y  '  tn  (^fc|fc  —  1  —  ■^■k\k—l)(^k\k—l  ~  l)  ) 

i 

where  Q  is  the  process  noise  covariance  matrix.  The  update 


The  variance  values  in  Q  associated  with  the  unknown  pa¬ 
rameters  determine  both  the  rate  of  parameter  estimation  con¬ 
vergence  and  the  estimation  performance  once  convergence 
is  achieved,  therefore,  techniques  have  been  developed  to 
tune  this  value  online  to  maximize  performance,  e.g.,  (Liu  & 
West,  2001;  Orchard,  Tobar,  &  Vachtsevanos,  2009;  Daigle, 
Saha,  &  Goebel,  2012).  We  adopt  the  approach  presented 
in  (Daigle,  Saha,  &  Goebel,  2012),  in  which  the  algorithm 
tries  to  control  the  variance  of  the  hidden  wear  parameter 
estimate  to  a  user-specified  range  by  modifying  the  process 
noise  variance.  Effectively,  the  algorithm  increases  the  vari¬ 
ance  when  the  relative  parameter  spread  is  below  the  desired 
level,  and  decreases  it  otherwise.  With  the  proper  settings, 
the  parameter  estimates  converge  quickly  and  track  with  high 
accuracy  and  precision. 

4.  Distributed  Prediction 

Each  local  prediction  module  takes  as  input  local  state- 
parameter  estimates  formed  from  the  local  estimators,  as  dis¬ 
cussed  in  Section  2.  The  required  estimates  must  be  con¬ 
structed  from  the  local  estimates  of  the  submodels  used  for 
estimation.  A  prediction  submodel  has  a  set  of  states  Xt 
and  parameters  ©j,  and  it  must  construct  a  local  distribu¬ 
tion  p(xj.,  9k\yo-k)-  To  do  this,  we  assume  that  the  local 
state-parameter  estimates  may  be  sufficiently  represented  by 
a  mean  /T  and  covariance  matrix  S\  For  each  prediction 
submodel  AT;,  we  combine  the  estimates  from  estimation 
submodels  that  estimate  states  and  parameters  in  Xt  U  B.( 
into  \xl  and  covariance  E\  If  there  is  overlap  in  the  state- 
parameter  estimates,  i.e.,  if  two  submodels  both  estimate  the 
same  state  variable  x  or  parameter  6,  then  this  may  be  re¬ 
solved  by  a  number  of  techniques,  e.g.,  taking  the  estimate 
with  the  smallest  variance,  or  taking  an  average.  Note  that, 
due  to  the  decomposition  into  independent  local  submodels. 
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Algorithm  1  EOL  Prediction 


Inputs:  {(x’(^,Vfc 


Outputs:  {EOLk 
for  j  —  I  to  .V  do 

k  <—  kp 

itt) 

A kP 
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V'--  )iwkP  Sj= 1 


*0)  .,AU)\N 


,w 


3  =  1 


‘(j) 
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°lU) 


Predict  uj. 
while  TlKOL(xk 
Predict  uj, 

Q‘U) 

& k.+  l 

Xfc+1  ~  P(x 
fe  -t—  fc  +  1 


O')  gi 0) 


,  ul)  =  0  do 


p(ei+iioi0)) 


Jfe+ll 


iOl^jO) 


0 


iO) 

K-fc 

,*0) 


end  while 


,*0) 

0*0) 

"fc+l 


EOL 

end  for 


*0) 


we  recover  only  an  approximation  to  the  joint  posterior  dis¬ 
tribution  as  would  have  been  found  by  a  global  estimator.  In 
particular,  covariance  information  is  lost  due  to  the  decou¬ 
pling  and  will  appear  as  zeros  in  the  merged  covariance  ma¬ 
trix.  As  shown  in  (Daigle  et  al.,  2011)  and  as  will  be  seen  in 
Section  5,  the  approximation  still  results  in  accurate  predic¬ 
tions. 

Given  the  mean  and  covariance  information,  we  represent  the 
distribution  with  a  set  of  sigma  points  derived  using  the  un¬ 
scented  transform.  Then,  as  in  (Daigle  &  Goebel,  2010),  each 
sigma  point  is  simulated  forward  to  EOL,  and  we  recover  the 
statistics  of  the  EOL  distribution  given  by  the  sigma  points. 

The  prediction  algorithm  is  executed  for  each  submodel  i,  de¬ 
riving  local  EOL  predictions  using  its  local  threshold  func¬ 
tion  based  on  the  local  EOL  constraints.  The  pseudocode  for 
the  prediction  procedure  is  given  as  Algorithm  1  (Daigle  & 
Goebel,  201  lb).  For  a  given  submodel  A4i,  each  sigma  point 
j  is  propagated  forward  until  T’EOL(x'k  '\  9'kJ> )  evaluates  to 
1.  The  algorithm  hypothesizes  future  inputs  uj. . 

Each  prediction  submodel  Al,  computes  a  local  EOL/RUL 
distribution,  i.e.,  p{EOL\p\yi0.kp)  and  p{RUL\p |yj:fcp). 
The  system  EOL  is  determined  by  the  minimum  of  all  the 
local  distributions,  since  Teol  of  the  system  is  1  whenever 
any  of  the  local  constraints  are  violated,  and  each  local  distri¬ 
bution  is  associated  with  a  subset  of  these  constraints.  Specif¬ 
ically,  for  to  prediction  submodels, 

p{EOLkp  |y0:fcp)  = 

To  compute  this,  we  sample  from  each  local  EOL  distribution 
and  take  the  minimum  of  the  local  samples.  This  is  repeated 
many  times  and  the  statistics  of  the  global  EOL  distribution 
are  computed. 


5.  Case  Study 

In  this  section,  we  apply  our  system-level  prognostics  ap¬ 
proach  to  a  four-wheeled  rover  testbed  developed  at  NASA 
Ames  Research  Center.  We  develop  a  model  of  the  rover,  and 
demonstrate  the  approach  using  simulated  scenarios. 

5.1.  Rover  Modeling 

The  rover  model  was  originally  presented  in  (Balaban  et  al., 
2011).  In  this  section  we  summarize  the  main  features  and 
include  some  extensions  to  the  model. 

The  rover  consists  of  a  symmetric  rigid  frame  with  four 
independently-driven  wheels.  The  wheel  speeds  are  governed 
by 

WFL  =  -7 ( TmFL  ~  TfFL  ~  TglFL  +  TgrFL )  (g) 

Jfl 

1 

WFR  =  7 -  (Tm.FR  —  TfFR  —  TgiFR  ~  TgrFRJ  (C2) 

■JFR 

WBL  =  7 ( TmBL  ~  TfBL  —  TgiBL  +  TgrBL )  (C3) 

Jbl 

UBR.  =  -J -  (TmBR  ~  TfBR  —  TgiFR.  ~  TgrBR )  (C4) 

•JBR 

The  F,  B,  L,  and  R  subscripts  stand  for  front,  left, 
back,  and  right,  respectively.  Here,  for  wheel  w  £ 
{FL,  FR,  BL,  BR },  ./„,  denotes  the  wheel  inertia;  Tmw  = 
kTiw  is  the  motor  torque,  where  iw  is  the  motor  current  and 
kT  is  an  energy  transformation  gain;  Tfw  =  p,fwojw  is  the 
wheel  friction  torque,  where  fi fw  is  a  friction  coefficient; 
Tgiw  =  rwpgis(vw  —  v)  is  the  torque  due  to  slippage,  where 
rw  is  the  wheel  radius,  pgis  is  a  friction  coefficient,  vw  is  the 
translational  wheel  velocity,  and  v  is  the  translation  velocity 
of  the  rover  body;  and  Tgrw  =  rU!ngrwui  cos  7  is  the  torque 
due  to  the  rotational  movement  of  the  rover  body,  where  pgrw 
is  a  friction  coefficient,  to  is  the  rotational  velocity  of  the  rover 
body,  and  7  =  arctan  l/b  with  l  being  the  rover  length  and  b 
being  its  width. 

We  consider  here  friction-based  damage  progression  in  the 
motors,  resulting  in  an  increase  in  motor  friction  over  time, 
which  will  lead  to  an  increase  in  power  consumption.  For 
wheel  w,  f.ifw  is  governed  by  (Daigle  &  Goebel,  201  lb) 


PfFL  =  VfFL  pfFL  W fl  (C5) 

PfFR  =  FfFR  PfFR  UpR  (cf) 

PfBL  =  FfBL  PfBL  U2bl  (c7) 

PfBR  =  VfBR  PfBR  ^br,  (C8) 


where  for  wheel  w,  Vfw  is  an  unknown  wear  coefficient. 

The  translational  velocity  v  of  the  rover  is  described  by 

v  =  —  ( FgiFL  +  FgiFR  +  FgiBL  +  ,  (C9) 
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Figure  2.  Rover  forces. 


where  m  is  the  rover  mass,  and  for  wheel  w,  Fgiw  = 
H giwiVw  ~  v)  is  the  force  due  to  slippage.  The  rotational 
velocity  oj  is  described  by 


The  temperature  of  the  windings  for  the  motors  are  governed 
by 


TdFL 

TdFR 

TdBL 

TdBR 


1 

Td 

1 

Td 

1 

Td 

l 

Td 


( i2FLR  —  hdFL{TdFL  —  TmFL )) 
(i2FRR  —  hdFR{TdFR.  —  TmFR)) 
{^BL^  hdBL(TdBL  ~  TmBL )) 
{}brR  ~  hdBR{TdBR.  —  TmBR ))  , 


(C19) 

(C20) 

(C2l) 

(C22) 


where  Jd  is  the  thermal  inertia  of  the  windings,  and  for  wheel 
w,  hdw  is  a  heat  transfer  coefficient,  and  Tmw  is  the  motor 
surface  temperature.  It  is  assumed  that  heat  is  lost  only  to 
the  motor  surface,  and  that  winding  resistance  R  is  approx¬ 
imately  constant  for  the  temperature  range  considered.  The 
surface  temperature  of  the  motor  for  wheel  w  is  given  by 


to  =  —(d  cos  7  FgiFR  +  d  cos  7  FgiBR  —  d  cos  7  FgiFL 

—  d  cos  7  FgiBL  —  dFgrFL  —  dFgrFR  —  dFgrBL 

—  dFgrBR).  (cio) 

Here,  J  is  the  rotational  inertia  of  the  rover,  d  is  the  distance 
from  the  center  of  the  rover  to  each  wheel,  and  for  wheel  w, 
Fgrw  =  AtgrwW  is  the  force  due  to  the  rotational  movement  of 
the  rover  body.  The  rover  forces  are  illustrated  in  Fig.  2. 

The  wheels  are  driven  by  DC  motors  with  PI  control  that  sets 
the  voltages  V  applied  to  the  motors.  The  motor  currents  are 
governed  by 


TmFL  =  ~r(hdFL(TdFL  ~  TmFL)  —  haFL{TmFL  ~  Ta)) 

J  S 

(C23) 

TmFR  =  ~r{hdFR{TdFR  —  TmFR )  —  haFR{TmFR  ~  Ta)) 

J  S 

(C24) 

TmBL  =  -y(hdBL{TdBL  —  TmBL )  —  haBL{TmBL  —  Ta)) 

j  s 

(C25) 

TmBR  =  -y{hdBR{TdBR  —  TmBR)  —  haBR{TmBR  —  Ta)), 

j  s 

(C26) 


iFL  =  ^TfL  —  iFLRFL  ~  k^U>FL)  (Cn) 

i-FR  —  Vfr  —  ^frRfr  —  kuujFR)  (ci2) 

=  YjTbL  —  IblRbL  —  kuUlBL)  (C13) 

IBR  =  -jr(VBR  —  IbrRbR  —  kuUBR)-  (C14) 


Here,  L  is  the  motor  inductance,  ku  is  an  energy  transforma¬ 
tion  term,  and  for  wheel  w,  R  is  the  motor  resistance.  The 
voltages  applied  to  the  motors  are  determined  by  the  con¬ 
trollers,  where  for  wheel  w,Vw  =  P  *  (uw  —  ojw)  +  I  *  e,iw, 
where  P  is  a  proportional  gain,  uw  is  the  commanded  wheel 
speed,  I  is  an  integral  gain,  and  is  the  integral  error  term. 


The  integral  error  terms  are  governed  by 

CiFL  =  UfL  —  UFL 

(Cl5) 

&iFR  =  UpR  ~  UpR 

(cie) 

eiBL  =  UBL  —  UBL 

(C17) 

CiBR.  =  UBR  ~  WBR- 

(ci8) 

The  motor  windings  heat  up  as  current  passes  through  them. 


where  Js  is  the  thermal  inertia  of  the  motor  surface,  and  for 
wheel  w,  haw  is  a  heat  transfer  coefficient,  and  Ta  is  the  am¬ 
bient  temperature.  Heat  is  transferred  from  the  windings  to 
the  surface  and  lost  to  the  environment. 

The  batteries,  which  are  connected  in  series,  are  described 
by  a  simple  electrical  circuit  equivalent  model  that  includes  a 
large  capacitance  Ci,  in  parallel  with  a  resistance  Rp,  together 
in  series  with  another  resistance  Rs  ,3  The  battery  charge  vari¬ 
ables  q,  are  governed  by 

qi  =  — Vi/Rpi  —  ( %fl  +  ifr  +  i'BR  +  ibl )  (C27) 

q_2  =  —  V2/R-P2  —  (iFL  +  iFR  +  i'BR  +  IBl)  (C28) 

93  =  —  Vz/RpZ  —  {iFL  +  IFR  +  IBR  +  IBl)  (C29) 

94  =  -  V4/Rp4  -  (iFL  +  iFR  +  iBR  +  Ibl)-  (c3q) 


3  We  use  a  simple  model  here  only  for  demonstration  purposes.  More  detailed 
battery  models  for  prognostics  can  be  found  in  the  literature,  e.g.,  (Saha  & 
Goebel,  2009). 
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M 10 

IFR,  eiFR 

0 

ufr,u*fr 

Cfr 

Cl2,Cl6,C37,C45 

Mu 

IBL,  e.iBL 

0 

UBL,U>bl 

Cbl 

Cl3,Cl7,C38,C46 

Mii 

%BR,  eiBR 

0 

UBR,R>BR 

Cbr 

Cl4,Cl8,C39,C47 

Mis 

LOFL,V,0J,  RfFL 

VfFL 

\flMfrMblMbr 

U*FL 

Cl,C5,C9,Cl0,C36,C45,C46,C47 

Ml4 

U>FR,V,U},  flfFR 

VfFR 

1fr,u}fl,u>blMbr 

UFR 

C2,C6,C9,Cl0,C37,C44,C46,C47 

M 15 

U>BL,V,U},  RfBL 

VfBL 

i blMflMfrMbr 

U*bL 

C3,C7,C9,Cl0,C38,C44,C45,C47 

M 16 

UlBR,  V,LU,RfBR 

VfBR 

i brMflMfrMbl 

Ubr 

C4,C8,C9,Cl0,C39,C44,C45,C46 

Table  1 .  Estimation  Submodels 


The  available  sensors  measure  the  voltages  of  the  batteries. 


given  as 


V-i  =  Qi/Cbi  —  Rsl  *  (*fl  +  *Ffi  +  +  *bl)  (C31) 

V2  =  q2/Cb2  —  RS2  *  (*FL  +  *Ffl  +  tBR  +  (C32) 

V3  =  q3/Cb3  —  Rs3  *  (*FL  +  *Ffl  +  IBR  +  * BL )  (C33) 

v;  =  qi/Cbi  -  f?s4  *  (*FL  +  *Ffi  +  IBR  +  IBl),  (C34) 

the  battery  current, 


ib  =  iFL  +  iFR  +  i'BR.  +  iBL, 

(C35) 

the  motor  currents. 

CpL  —  iFL 

(C36) 

i*FR  =  iFR 

(C37) 

i*BL  =  iBL 

(C38) 

i*BR  =  iBR, 

(C39) 

the  motor  surface  temperatures. 

/Ti*  rji 

1 mFL — 1mFL 

(C40) 

FTl*  rji 

1 mFR  —  1 mFR 

(C41) 

ft~i*  rji 

1 mBL  —  4  mBL 

(C42) 

rj~i*  rji 

1  mBR  —  1  rn.BR, 

(C43) 

and  the  wheel  speeds. 

WFL  =  WFL 

(c44) 

U>FR  =  wfr 

(C45) 

UBL = ^BL 

(C46) 

UBR  ~  UBR- 

(C47) 

Ei  >  V~ 
E2  >  V~ 

E3  >  v- 

E4  >  V~ 


(c48) 

(C49) 

(C50) 

(C51) 


where  the  voltage  threshold  is  given  by  V  -  9.6  V,  and  for 
battery  i,  V)  =  qi/Cbi-Rsi*(iFL+iFR+iBR+iBL)-  We  are 
also  interested  in  when  the  motor  temperature  gets  too  high. 
The  motor  windings  are  designed  to  withstand  temperatures 
up  to  a  certain  point,  after  which,  the  insulation  breaks  down, 
the  windings  short,  and  the  motor  fails  (Balaban  et  al.,  2010). 
The  constraints  are  given  as 


TmFL  <  T+ 
TmFR  <  T m 
TmBL  <  T+ 
TmBR  <  T+ , 


(C52) 

(c53) 

(C54) 

(C55) 


where  the  temperature  limit  is  given  by  7j+  =  70°  C.  The 
rover  cannot  be  operated  when  any  of  these  constraints,  c48~ 
C55,  are  violated. 

In  the  general  case,  we  consider  uncertainty  in  the  friction 
wear  parameters  VfpL ,  vfFR,  f/bl,  and  f/br',  the  heat 
transfer  coefficients  hdFL »  hdFR ,  hdBL,  hdBR »  hapL »  haFR, 
haBL ,  and  haBR :  the  battery  capacitances  Cbi,  Cb2 ,  Cft3 ,  and 
Cb4\  and  the  battery  resistances  Rs  1,  f?s2,  Rs3,  and  f?s4.  Sen¬ 
sor  and  process  noise  were  estimated  based  on  data  from  the 
actual  rover  testbed. 


Here,  the  *  superscript  indicates  a  measured  value. 

We  are  interested  in  predicting  when  any  of  the  rover  batter¬ 
ies  are  at  their  voltage  threshold,  beyond  which  the  batteries 
will  be  damaged  (Saha  &  Goebel,  2009).  The  constraints  are 


5.2.  Results 

To  demonstrate  the  validity  of  the  approach,  we  describe 
two  scenarios  for  system-level  prognostics  of  the  rover.  In 
the  first,  the  rover  is  operating  nominally  without  any  faults 
present,  and  in  the  second,  friction  damage  is  progressing 
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Submodel 

Xl 

<5/ 

L h 

Yi 

AI17 

Qi 

Cbl,  Ksi 

lFL,iFR,  i-BL,  IBR 

0 

C27.C48 

M 18 

<72 

Cb2 ,  Rs  2 

IFlUfR,  i-BL,  IBR 

0 

C28.C49 

M19 

<73 

Cb3,  Rs3 

lFL,iFR,  tBL ,  IBR 

0 

C29.C50 

M20 

94 

f)>4  >  Rs4 

IflUfr,  iBL,  iBR 

0 

C3O.C51 

M21 

TdFL,TmFL 

hdFL,  haFL 

IFL 

0 

Cl9.C23.C52 

M22 

TdFR>  TmF  R 

fodFRi  ha.FR 

IFR 

0 

C20.C24.C53 

M23 

TdBL,TmBL 

hdBL,  haBL 

i-BL 

0 

C21,C25.C54 

M24 

LdBRi  TmBR 

hdBR,  haBR 

IBR 

0 

C22.C26.C55 

Table  2.  Prediction  Submodels  Using  Motor  Currents  as  Local  Inputs 


on  one  motor.  In  both  cases,  the  rover  travels  between  var¬ 
ious  waypoints,  moving  at  an  average  speed  of  0.5  m/s.  The 
unknown  parameters  are  initialized  incorrectly  (with  around 
10%  error)  so  that  the  local  estimators  must  converge  to  the 
true  values.  In  both  cases,  the  estimation  step  is  performed 
in  a  distributed  manner  using  the  set  of  submodels  derived 
by  using  measured  values  as  local  inputs,  shown  in  Table  1 . 
For  example,  submodel  A4 1  computes  an  estimate  of  V*  us¬ 
ing  the  measured  value  of  fj*  as  a  local  input,  and  using  the 
minimal  set  of  constraints  to  do  this.  For  the  prediction  sub¬ 
models,  as  will  be  shown,  the  correct  submodels  to  use  de¬ 
pends  on  the  scenario,  and  illustrates  when  and  when  not  the 
prediction  step  can  be  decomposed. 

5.2.1.  Nominal  Operation 

We  first  consider  a  scenario  involving  nominal,  fault-free  op¬ 
erations.  In  this  case,  EOL  occurs  around  3  h.  An  RUL  pre¬ 
diction  is  made  every  500  s.  With  the  rover  traveling  at  an 
average  speed  of  0.5  m/s,  the  motor  currents  average  to  about 
0.15  A  each  and  so  the  total  current  draining  from  the  four 
batteries  is  0.6  A.  Since  these  values  do  not  vary  much  dur¬ 
ing  nominal  operation,  we  can  use  the  motor  currents  as  local 
inputs  for  the  model  decomposition.  These  submodels  are 
shown  in  Table  2.  Note  that  the  estimates  from  the  estima¬ 
tion  submodels  A4i-A4g  are  used  directly  in  the  prediction 
submodels  A4 17- A4 24,  respectively,  and  that  estimation  sub¬ 
models  AI9-AI16  are  not  necessary.  Note  also  that  the  pre¬ 
diction  submodels  do  not  compute  any  outputs,  rather,  their 
goal  is  to  compute  EOL  constraints  (e.g.,  A4  \  7  computes  c,^). 

The  system-level  prediction  results  are  shown  in  Fig.  3.  Pre¬ 
dictions  from  the  battery  submodels  are  shown  in  Fig.  4.  In 
this  case,  the  motor  temperatures  reach  a  steady-state  that  is 
below  the  temperature  threshold,  so  only  the  batteries  have  an 
impact  on  system  EOL,  which  is  the  minimum  of  the  EOLs 
predicted  for  the  battery  submodels.  In  particular,  it  is  the 
first  and  fourth  batteries  (corresponding  to  AI17  and  .VI 20, 
respectively)  that  discharge  the  fastest,  as  shown  explicitly  in 
Fig.  4.  The  figures  show  the  means  of  the  predicted  RUL 
distributions,  the  true  RUL,  RUL* ,  and  an  accuracy  cone  of 
a  =  10%  around  it.  In  Fig.  3,  we  show  both  the  system- 
level  predictions  using  the  distributed  approach  with  A4 17- 
Ad  24  and  the  centralized  approach  using  the  global  predic¬ 
tion  model  A4q-  The  global  prediction  model  contains  all 


the  states,  parameters,  and  constraints  given  in  the  previous 
subsection,  minus  the  output  constraints,  and  uses  the  com¬ 
manded  wheel  speeds  (known  a  priori)  as  hypothesized  in¬ 
puts.  Since  the  currents  are  also  known  a  priori,  the  system- 
level  prediction  can  be  decomposed,  and  the  predictions  made 
using  the  local  submodels  closely  match  those  made  using  the 
global  model,  as  shown  in  the  figure. 

We  use  the  relative  accuracy  (RA)  metric  (Saxena  et  al.,  2010) 
for  prediction  accuracy.  Averaged  over  all  predictions,  RA  is 
97.48%  for  the  distributed  approach  and  98.74%  for  the  cen¬ 
tralized  approach.  Using  relative  standard  deviation  (RSD)  as 
a  measure  of  spread,  and  averaged  over  all  prediction  points, 
RSD  is  0.40%  for  the  distributed  approach  and  0.43%  for 
the  centralized  approach.  The  distributed  approach  is  only 
slightly  less  accurate  but  has  better  precision.  Here,  both  ap¬ 
proaches  are  very  accurate  since  the  system  state-parameter 
estimates  are  very  accurate,  and  there  is  only  a  small  amount 
of  error  associated  with  assuming  a  constant  average  mo¬ 
tor  current  or  wheel  speed.  Correspondingly,  the  prediction 
spread  is  relatively  small  because  the  uncertainty  in  the  state- 
parameter  estimate  is  very  small. 

5.2.2.  Friction  Damage  Progression 

We  now  consider  a  scenario  in  which  for  the  front-left  motor, 
there  is  nonlinear  friction  damage  progression  with  vjfl  = 
1  x  10-4  s.  As  a  result  of  the  continuously  increasing  fric¬ 
tion,  the  current  drawn  by  the  motor  increases  as  well  in  order 
for  the  motor  controller  to  maintain  the  same  desired  wheel 
speed.  Hence,  the  total  current  drawn  from  the  batteries  is 
increased,  and  EOL  occurs  around  2  h.  Because  ip l  is  con¬ 
stantly  changing,  and  in  a  way  that  is  dependent  on  the  motor 
state,  it  cannot  be  predicted  a  priori ,  and  so  cannot  be  used  as 
a  local  input  because  the  resulting  predictions  will  not  be  ac¬ 
curate.  Therefore,  we  require  a  submodel  that  estimates  iph , 
and  we  so  employ  submodels  using  as  local  inputs  average 
values  for  the  remaining  motor  currents,  average  commanded 
wheel  speeds,  and  average  rover  translational  velocity  v  and 
rotational  velocity  to.  The  prediction  submodels  for  this  case 
are  shown  in  Table  3.  For  comparison,  we  demonstrate  also 
prediction  using  A4i7-A424,  and,  for  this  strategy,  at  each 
prediction  point  the  average  value  of  current  measured  over 
the  last  minute  is  used  as  the  future  hypothesized  value.  Of 
course,  this  will  not  yield  accurate  results  since  future  values 
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Submodel 

Xi 

65) 

Ui 

Y~ 

M25 

qi,lFL,  eiFL ,  U>FL ,  HfFL 

U'fti,  Hsl,  F  fpL 

UFL,V ,  L 0,  IFR,  lBL,lBR 

0 

Cl,C5,Cll,Cl5,C27,C48 

Ad  26 

q2,lFL,  eiFL ,  U>FL  ,  fJfFL 

Cb2,  Rs2,  FfFL 

UFL,V,  UJ,  IFR,  iBL,lBR 

0 

Cl ,C5,Cl 1 ,Cl5 ,C28 ,C49 

AA.27 

?3 ,  * FL ,  eiFL ,  UJFL  MfFL 

Ub3,  Rs3,  FfFL 

Cb4,  Rs 4,  FfFL 

UFL,V,  UJ,  IFR,  iBL,lBR 

0 

Cl  ,C5,Cl  1  ,Cl5  ,C29,C50 

M2g 

q4,,lFL,  eiFL ,  UFL  MfFL 

UFL,V,  UJ,  IFR,  iBL,lBR 

0 

Cl,C5,Cn,Cl5,C30,C5l 

Ad  29 

TdFL ,  TmFL  ,lFL,eiFL,  U >FL ,  ft fFL 

hdFL,  haFL ,  FfFL 

UFL,V,UJ 

0 

Cl9,C23,C52,Cn,Cl5,Ci,C5 

Ad  30 

TdFRi  TmFR 

hdFRi haFR 

IFR 

0 

C20,C24,C53 

Ad3l 

TdBL,TmBL 

hdBL,  haBL 

iBL 

0 

C21,C25,C54 

Ad  32 

TdBR ■>  TmBR 

hdBR, haBR 

iBR 

0 

C22,C26,C55 

Table  3.  Prediction  Submodels  Using  Commanded  Wheel  Speeds  and  Rover  Velocities  as  Local  Inputs 
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Figure  3.  System  RUL  prediction  performance  under  nominal 
conditions  with  a  =  0.1. 


Figure  4.  Individual  battery  submodel  RUL  prediction  per¬ 
formance  under  nominal  conditions  with  a  =  0.1. 


of  the  current  will  actually  be  larger.  Note  that  the  predic¬ 
tion  submodels  used  in  this  case  do  not  correspond  directly  to 
those  used  for  estimation.  So,  when  constructing  the  estimate 
for  A/1 25,  for  example,  it  takes  the  estimates  from  Adi,  Ad 9, 
and  Ad  13. 

The  system-level  prediction  results  are  shown  in  Fig.  5.  Al¬ 
though  the  increased  friction  causes  the  temperature  of  the 
front-left  motor  to  increase,  it  is  still  the  batteries  discharging 
that  dominates  the  system-level  EOL  in  this  case.  We  show 
the  predictions  using  Adi7-Ad24,  M.23-M-32,  and  the  global 
model  Ado-  For  M25-M32,  average  values  of  v  =  0.5  m/s 
and  u>  =  0  rad/s  are  used.  Here,  the  predictions  using  the  lat¬ 
ter  two  approaches  are  virtually  identical  (the  predictions  us¬ 
ing  Ad  25-Ad  32  are  hidden  under  those  for  Ado),  and  fairly  ac¬ 
curate.  In  contrast,  as  expected,  the  predictions  using  Ad  17- 
Ad  24  are  very  inaccurate,  and  only  converge  towards  the  true 
RUL  at  the  very  end.  This  quite  effectively  demonstrates 
that,  in  this  scenario,  it  is  incorrect  to  use  the  front-left  mo¬ 
tor  current  as  a  local  input  for  predictions,  since  it  cannot  be 
predicted  independently  from  the  front-left  motor  submodel. 


and  therefore  a  submodel  that  itself  predicts  this  current  is 
required  to  obtain  accurate  predictions. 

Here,  RA  averages  to  58.95%  using  Ad  17-Ad 24,  94.24%  us¬ 
ing  Ad25-Ad32,  and  94.32%  using  Ado-  RSD  averages  to 
0.76%  using  Adi7-Ad24,  1-62%  using  Ad25-Ad32,  1.73% 
using  Ado-  Here,  we  also  observe  an  increase  in  prediction 
spread  using  the  centralized  approach  with  only  a  slight  in¬ 
crease  in  accuracy  over  the  distributed  approach.4  Overall, 
accuracy  and  precision  are  both  decreased  compared  to  the 
nominal  scenario  because  there  is  more  uncertainty  in  the 
state-parameter  estimate,  specifically,  that  dealing  with  the 
estimate  of  z//fl-  This  uncertainty  in  the  state-parameter  es¬ 
timate  contributes  to  the  additional  uncertainty  in  the  RUL 
predictions. 


4The  RSD  for  A4i7-A424  is  the  lowest  because  those  submodels  do  not 
include  the  motor  friction  component,  so  do  not  have  the  additional  uncer¬ 
tainty  associated  with  the  estimation  of  the  wear  par  ameter. 
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Figure  5.  System  RUL  prediction  performance  with  friction 
damage  progression  with  a  =  0.1. 


6.  Conclusions 

In  this  paper,  we  formulated  the  system-level  prognostics 
problem  and  proposed  a  distributed  solution  based  on  struc¬ 
tural  model  decomposition.  Using  a  four-wheeled  rover  as 
a  simulation-based  case  study,  we  demonstrated  the  effec¬ 
tiveness  of  the  approach.  Most  importantly,  the  distributed 
approach  allows  for,  in  many  practical  circumstances,  the 
decomposition  of  the  system-level  prognostics  problem  into 
component-level  prognostics  problems  and  provides  a  mech¬ 
anism  to  merge  local  prognostics  results  into  a  system-level 
result.  Further,  since  the  local  subproblems  are  independent, 
this  allows  component  experts  to  focus  on  prognostics  solu¬ 
tions  for  their  components.  However,  we  showed  also  that 
this  approach  is  not  always  possible  if  accurate  results  are  de¬ 
sired,  since  in  some  cases  the  prediction  problem  cannot  be 
so  easily  decomposed,  and  it  depends  crucially  on  correct  as¬ 
sumptions  about  what  variables  may  serve  as  local  inputs  for 
the  prediction  problem. 

Although  in  this  paper  we  focused  on  the  model-based  prog¬ 
nostics  paradigm,  our  approach  is  flexible  in  that  data-driven 
algorithms  may  be  used  also,  once  the  local  subproblems  are 
defined.  For  example,  in  previous  works,  structural  model 
decomposition  was  used  to  automatically  design  gray  box  di¬ 
agnosis  models  that  were  implemented  using  different  data- 
driven  techniques  (for  instance,  state  space  neural  networks 
in  (Pulido,  Zamarreno,  Merino,  &  Bregon,  2012)  or  machine 
learning  techniques  in  (Alonso-Gonzalez,  Rodriguez,  Prieto, 
&  Pulido,  2008)).  By  decomposing  the  system-level  prob¬ 
lem  into  independent  subproblems  through  structural  model 
decomposition,  we  can  apply  similar  ideas  to  solve  each  prog¬ 


nostics  subproblem  by  using  the  most  appropriate  technique, 
whether  it  is  a  model-based,  data-driven,  or  hybrid  approach. 

An  important  direction  of  future  work  is  in  algorithms  for 
optimal  placement  of  sensors  for  model  decomposition,  be¬ 
cause  the  level  of  model  decomposition  that  can  be  achieved 
is  dependent  on  the  number  of  sensors  and  where  they  are 
placed.  This  results  in  the  most  efficient  decomposition  of 
the  system-level  prognostics  problem.  Current  work  also 
addresses  combining  the  distributed  prognostics  framework 
with  a  distributed  diagnostic  approach  for  integrated  diag¬ 
nostics  and  prognostics  (Bregon,  Daigle,  &  Roychoudhury, 
2012). 
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